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ABSTRACT 



The emergence of magnetic flux from the convection zone into the corona is an important 
process for the dynamical evolution of the coronal magnetic field. In this paper we extend our 
previous numerical investigations, by looking at the process of flux interaction as an initially 
twisted flux tube emerges into a plane parallel, coronal magnetic field. Significant differences 
are found in the dynamical appearance and evolution of the emergence process depending on 
the relative orientation between the rising flux system and any preexisting coronal field. When 
the flux systems are nearly anti-parallel, the experiments show substantial reconnection and 
demonstrate clear signatures of a high temperature plasma located in the high velocity outflow 
regions extending from the reconnection region. However, the cases that have a more parallel 
orientation of the flux systems show very limited reconnection and none of the associated features. 
Despite the very different amount of reconnection between the two flux systems, it is found that 
the emerging flux that is still connected to the original tube, reaches the same height as a function 
of time. As a compensation for the loss of tube flux, a clear difference is found in the extent of the 
emerging loop in the direction perpendicular to the main axis of the initial flux tube. Increasing 
amounts of magnetic reconnection decrease the volume, which confines the remaining tube flux. 

Subject headings: Sun: magnetic fields - Numerical experiments - Sun: active regions - Sun: corona 



1. Introduction 

Flux emergence is one of the manifestations 
of the continuously changing solar magnetic field. 
In this process magnetic flux is transported up 
through the convection zone presumably by a 
combination of buoyancy and convection towards 
the photosphere. At the photosphere the phys- 
ical structure of the sun changes, from a con- 
vectively unstable to a convectively stable atmo- 
sphere, making the continued rise of magnetic flux, 
due to buoyancy, more difficult. Despite this, nu- 
merous observations in all wavelength ranges show 
that the emergence process is a frequently occur- 



ring event in the solar atmosphere. 

There are numerous observations of the pho- 
tosphere and lower transition region (iLites et al 



1995 



1998 : IStrouse fc Zwaanl [l999; Kubo et al. 



20031 : IPariat et al. 2004; Lites 2005). These show 
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picture whereby emerging magnetic fields change 
the convective flows, allowing for initially horizon- 
tal magnetic fields to penetrate the photosphere. 
As the magnetic field expands into the transition 
region and lower corona, relatively cold plasma is 
lifted by the magnetic field and eventually starts 
draining back towards the photosphere, along the 
magnetic field lines. In isolated emergence events 
two strong opposite polarities arise forming two 
magnetically connected sunspots. In the region 
between the two flux concentrations new flux con- 
tinues to emerge, with new positive - negative 



pairs arising at different positions. iLites et al 
(|1995l ) used observations covering temperatures 
from photosphere to corona plasma to argue that 
a full coherent flux tube rises from the convection 
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zone. IStrouse fc Zwaanl ( 1999f ) suggested that the 
flux emerges in "parallel" sheets between the two 
major flux concentrations, and that the random 
appearance of new small flux concentrations be- 
tween the two sunspots favor a picture whereby 
undulating field lines emerge at different locations 
along the full leng th of the fiel d line connecting 
the two sunspots. IPariat et al.l (|2004l ) uses field 
extrapolation of the photospheric field to estab- 
lish a field line structure of the magnetic field and 
find that undulating field lines exist in this model. 
Furthermore , the dips are associat ed with Eller- 
man bombs ( Georgoulis et al.l 12002). and they in- 
terpret these as being responsible for leaving dense 
material in the part of the loop that dips below 
the photosphere, making the emergence for the re- 
maining upper part of the loop much easier. 

The manifestation of flux emergence in the so- 
lar corona often results in a strong interaction be- 
tween the two initially disconnected flux systems. 
This interaction results in local brightenings ob- 
served in all wavelengths from white light through 
EUV and into X-ray. Often these events are as- 
sociated with high-speed flows of the hot plasma 
emerging from the (reconn ection) region where the 
two flux systems interact. lLongcope et al.l (j2QQ5h 
investigated one such event. Using simple flux es- 
timates for the flaring structures, they showed how 
the new flux regions changed its connectivity with 
time and became connected to the existing coronal 
magnetic field. This happened in such a way that 
no dynamic activity took place for a initial period 
of time, followed by a rapid change of connectivity, 
and finally a more quiet phase. 

How does this interaction depend on the struc- 
ture of the two initially separate flux systems? No 
clear analysis has been made on existing observa- 
tional data, and one can ask if this is actually pos- 
sible, with the present inability of directly tracking 
magnetic field lines and their photospheric connec- 
tions. Further to this, each observed event is dif- 
ferent from previous ones in a number of ways, so 
how can one quantify the most important reasons 
for these differences in evolution? 

Such effects are much easier to investigate by 
undertaking a series of numerical experiments 
where one can control, in detail, the environment 
into which t he flux emergence occurs. In previous 
experiments ([F an 2001; Ma gara fc Longcope1l2003l: 
Archontis et al.1 I2QQ4L 120051 : iManchester et al.1 



2004 : Magara 20061 ). the rise of a twisted loop is 
initiated inside the convection zone and followed 
in a self-consistent manner. In these investiga- 
tions, the structure of the emerging flux resembles 
a situation where the outer layers of the flux tube 
expand into the corona. In a manner, this resem- 
bles the process of peeling off the outer layers of an 
onion and e xpanding these parts into the corona. 
Specifically iMagaraT ( 20061 ) . compared the devel- 
opment of the photospheric flux concentrations 
with observations, finding the same characteris- 
tic evolution of the flux concentrations and their 
associate d inversion polar i ty line as seen in obser- 
vation s. iGalsgaard et al.l (l2QQ5h : I Archontis et al 



(|2005l ) and llsobe et al.l (|2005l ) presented the first 
3D MHD experiments in which a twisted flux tube 
emerges from the convection zone and interacts 
with a pre-existing coronal magnetic field. 



Archontis et al.l (|2QQ5h analyzed in detail the 



interaction between an emerging magnetic flux 
system and a uniform horizontal coronal magnetic 
field, using for the rising tube the same i nitial 
condit i ons as in previous exp eri ments (e.g. iFan 



200lh : lArchontis et al.1 (|20Q4h ). lArchontis eTal 



20051 ) show how the sub-photospheric flux tube 



emerges into the corona and pushes the magnetic 
field upward and outward. Given the initial al- 
most antiparallel, mutual orientation of the sys- 
tem at the time of first contact, a strong current 
sheet is formed at the interface. The interaction 
of the two flux systems then follows a complicated 
pattern that slowly changes in time. This inter- 
action depends also on the relative orientation be- 
tween the two sy s tems, as has been pointed out by 
Galsgaard etall (|2QQ5h . 



In a different approach, IFan fc Gibson! (12004) 
used a 2D magnetic arcade field embedded in a 
constant temperature coronal region, into which 
they force a pre-defined curved loop structure 
using an imposed boundary electric field. This 
showed, for the setup with minimal reconnection 
between the two flux systems, that the twisted 
emerging loop entered into the "corona", even- 
tually experienced a kink instability and strong 
currents were generated where reconnection takes 
place due to non-ideal effects. The emergence 
of a twisted l o op in to the corona is supported 
by iLites et aL (1995). Therefore it is interest- 
ing to compare the structure of curr ent concen- 
trations between [Fan fc Gibson! ( 20041 ) and exper- 
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iments which included the full emergence process, 
knowing that differences in the coronal structure 
and the emergence process may provide very dif- 
ferent dynamical evolutions. 

In thi s paper, we follow up on t he work pre- 
sented in lGalsgaard et al] (|2005f ) and lArchontis et al 
(|2QQ5l ) and investigate, in more detail, the signif- 
icant deviations in the emergence process as the 
orientation between the two interacting magnetic 
flux systems is changed. In Section [2] we briefly 
describe the numerical model. Sections 3-6 con- 
tain the various results of the experiments and 
Section [71 is a discussion of the implications of 
these results. Finally, highlights of the results are 
summarized in Section El 

2. Model setup and numerical approach 

The parameters of the magnetic flux tube, the 
background stratification and t he initial condi- 



tions f ollo w the work presented inlGalsgaard et al 



(|2QQ5l ) and lArchontis et al.1 (|2QQ5l ). As a reference 
for the following discussion we start by summa- 
rizing the initial state of the experiments and the 
numerical approach. 

Our model consists of a highly stratified en- 
vironment and a horizontally twisted magnetic 
tube. The background medium consists of an adia- 
batically stratified convection zone, an isothermal 
layer representing the photosphere, a region where 
the temperature steeply increases with height and 
represents the transition region and finally an 
isothermal layer with coronal temperatures. The 
tube center is located almost 2 Mm below the base 
of the photosphere. The longitudinal component 
of the magnetic tube has a Gaussian profile with a 
central field strength of 3.8 kG, while the twist is 
uniform around the axis of the tube. This particu- 
lar flux tube is stable towards the kink instability. 
This gives a plasma (3 = 12.8 at the axis of the 
tube. The rise of the tube is triggered by imple- 
menting a density deficit distribution that has a 
maximum at the middle of the axis of the tube 
and with a Gaussian distribution along the tube. 
Initially, a horizontal magnetic field is included in 
the atmosphere above the lower transition region. 
The orientation of the ambient field relative to the 
main axis of the tube is an important parameter 
in the experiments. 

Figure [T] shows the gas pressure, temperature 
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Fig. 1. — Distribution of gas pressure (thick 
solid), density (dashed), temperature (dash- 
dotted) and magnetic pressure (thin solid) along 
the central, vertical (y = 0,x = 0) line. 



and density of the stratified environment as a func- 
tion of height. All the profiles are normalized ac- 
cording to the photospheric values: p h = 1.4 10 5 



erg cm 



-3. 



Pph 



= 310- 7 



g cm 



" 3 ; T h = 5.6 10 3 



K and H h = 170 km. Other units used in the 
simulations are: time, t h = 25 sec; velocity, 
V = (p p h/ Pph) 1 ^ 2 = 6-8 km sec -1 and magnetic 
field, B ph = 1.3 10 3 Gauss. 

The distribution of the magnetic pressure in 
Fig. [1] shows the magnetic flux tube and the am- 
bient field. The ambient field is given by 



B c (z) [cos(0),sin(0),O], 



(1) 



where B c (z) is described by an hyperbolic tangent 
profile. The intensity of the coronal field is chosen 
such that the local plasma (3 is close to 0.06. 

The direction of the initial ambient field is given 
by the polar angle (/>, which is measured in a hor- 
izontal x?/-plane from the positive x— axis. At the 
same time the magnetic fleldlines at the top of 
the rising flux system are oriented in an approxi- 
mately antiparallel direction to the ambient field- 
lines when (j) = and they are almost parallel 
when (j) = 180 (see Fig.[2j). Thus, the relative hor- 
izontal angle between the two flux systems upon 
contact is 0o = 180 — <p deg. The polar angle, 0, 
of the coronal magnetic field is different in the five 
experiments we have performed; it changes from 
(j) = to (j) = 180 (see Table (pQ)). 
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Fig. 2. — Top view of an emerging twisted fieldline and the ambient coronal field (long horizontal arrows). 
The short arrow at the center of the box shows the direction of the emerging fieldline. The three panels 
represent different relative orientations of the two flux systems: experiment A: (j) = 0, experiment B: (j) = 45 
and experiment C: <\> = 90. 



Experiment 


00 


Line style 


Orientation 


A 


180. 


Triple Dotted dashed 


Antiparallel 


B 


135. 


Dot-Dashed 


Slanted 


C 


90. 


Full 


Perpendicular 


D 


45. 


Long-Dashed 


Slanted 


E 


0. 


Dotted 


Parallel 



Table 1: List of the numerical experiments performed in the simulations. The first column shows the reference 
names of the experiments. The second column shows the relative horizontal angle of the two magnetic 
systems. The third column indicates the linestyle used in the different plots in this paper and the fourth 
column the orientation of the coronal field relative to the upcoming field. 



The evolution of the system is governed by 
the three-dimensional, time-dependent and resis- 
tive MHD equations. These are solved using a 
numerical approach based on high order finite dif- 
ferencing on staggered grids. By using 6 neigh- 
boring data points, a 6th order accurate spatial 
derivatives and corresponding 5th order accurate 
interpolations routines are used. The solution 
is advanced in time using a 3rd order predictor- 
correction algorithm. Due to the high spatial or- 
der, special treatment of viscosity and resistivity 
are required to prevent numerical ringing in the 
vicinity of steep gradients in the physical quanti- 
ties. This is handled by a combined approach, that 
is designed to remove numerical problems that can 
occur in specific problems. These approaches are 
localized in space, implying that dissipation only 
take place over length scales of a few gridpoints. 
Using such an approach, makes it impossible to 



assign a single c haracteristic Reynolds numbe r for 
the experiment ( Nordlund fc Galsgaardl 1997 ). 



The numerical resolution of the experiments is 
(148, 160, 218) in the (x, y, z) directions, with z be- 
ing the height. The size of the numerical domain is 
(-60, 60), (-70, 70) and (-22, 70), which is equiv- 
alent to a box of sides 20.4 Mm x 23.8 Mm x 15.6 
Mm. The resolution in the x and y directions is 
137.8 km/cell and 148.8 km/cell correspondingly. 
The grid in the vertical direction is stretched in 
a way that the highest resolution covers the re- 
gion from the top of convection zone to the bot- 
tom of the corona. Here the grid resolution is 47.7 
km/cell. The resolution has lower values close to 
the top and bottom boundary of the numerical 
domain. 
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3. Current sheet orientation 

The location, orientation and strength of cur- 
rent sheets in 3D is vital for providing an environ- 
ment for fast magnetic energy release. This sec- 
tion is concerned with investigating the buildup 
of current sheets in the various experiments, with 
the aim to study the relation between some basic 
model parameters and locations of reconnection. 

As the emerging flux pushes its way into the 
coronal magnetic field, stress builds up at the in- 
terface between the two flux systems. When the 
two flux systems are antiparallel (<fo = 180), a cur- 
rent concentration is formed all over the emerging 
plasma hill. As the stress continues to build up the 
current is concentrated into a narrow curved sheet 
that reaches from the summit point of the plasma 
hill down its sides towards the photosphere al- 
most along the direction of the underlying emerg- 
ing flux tube in the ^-direction. As the orientation 
of the coronal magnetic field changes in the dif- 
ferent experiments, the orientation and strength 
of the current sheet changes too. The current 
sheet is found to rotate around its vertical cen- 
tral axis as a monotonic function of the angle 
</>o between the coronal magnetic field and the 
emerging magnetic field (see also Fig. [3]), where 
cos 0o = (B cor • Bt)/|B cor ||Bt|, B t represents the 
tube field at the summit point of the emergence 
region. 



<fio (degrees) 

Fig. 3. — The orientation of the current sheet with 
respect to the vertical yz midplane as a function of 
the angle between the two magnetic flux systems. 

The orientation of the current sheet can be 
found using the following analysis. Assume we are 
only interested in the orientation at the summit 



point of the emergence, then the z component of 
the field can be ignored. In general, the field vec- 
tors can be expressed in terms of an orthogonal co- 
ordinate system, with one unit vector, ei, defined 
by the direction given by the sum of the magnetic 
vectors, B cor + B^ and the other unit vector, e2, 
orthogonal to this defining a right hand system. 
Thus, 



B t 



B cor \\ e i + B cor ±e2, (2) 
Bt||ei + Bt±e 2 (3) 



The components of B cor and B t along ei, namely 
B cor \\ and B t \\, represent the components of the 
magnetic field that cannot be annihilated and at 
the same time determine the direction of the main 
axis of the current sheet. The components along 
^2, B cor ±_ and B t ±, provide two oppositely di- 
rected components that can annihilate in a recon- 
nection process. The integrated current across the 
sheet is simply given by \B cor ± \ + |-Bt_i_|. Visual- 
ization of the current sheet and illustration of its 
orientation in the computational volume is shown 
in Section [6l 



3.1. Magnetic Pressure 
the Current Sheet 



Balance Across 



In 2D the total pressure balance across a cur- 
rent sheet is simple, and it maintains a change 
between magnetic and gas pressure. In 3D such a 
balance may change with time as the relative ori- 
entation of the field that is advected into the sheet 
changes. In this paper, the evolution of the total 
pressure balance is used as an indirect indicator of 
magnetic reconnection. 

As the tube rises the current sheet moves up- 
ward co vering the upper part o f the buoyant flux 



system. lArchontis et al.l (|2QQ5[ ) showed that the 



magnetic pressure distribution across the current 
sheet changes with time in experiment A. This 
change in the temporal evolution of the magnetic 
pressure was found to play an important role in 
the reconnection process between the two flux sys- 
tems. 

Figure [4] shows the pressure of the horizon- 
tal component of the magnetic field (Bh or ) below 
(solid line, |B t |) and above the current sheet (dot- 
dashed line, |B cor |) as a function of time. The 
dashed line represents the pressure from the B x 
component of the magnetic field below the current 
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Fig. 4. — Temporal evolution of the magnetic pressure just above and below the current sheet are shown 
for experiments A, B, C and D from top- left to bottom-right panel. Full line represents the pressure below 
the sheet, the dot dashed line indicates the pressure above the sheet. Finally the dashed line indicate the 
pressure from the magnetic component perpendicular to the tube direction just below the sheet. 



sheet. 

At the beginning of the emergence process the 
magnetic pressure inside the rising plasma is much 
larger than that of the ambient coronal field. At 
t ~ 60 the difference is more than one order 
of magnitude. This pressure excess pushes the 
emerging tube upwards into the atmosphere - for 
details see the discussion of forces given in Sec- 
tion I5.21 After t « 80 a transverse balance of 
magnetic pressure is achieved, independently of 
the different orientation of the ambient field, and 
this balance remains until the end of the simula- 
tion. 

Before t ~ 95 the B x and the Bhor of the ris- 
ing magnetic field in the tube, are approximately 
equal. For experiments A and B the two compo- 
nents separate around this time, indicating that 
the orientation of the emerging magnetic field just 



below the interface changes with time. The reason 
being that, as time proceeds, the uppermost field 
lines of the rising plasma reconnect with the ambi- 
ent field and this allows for different internal flux 
layers to come into contact with the overlying flux 
system; the magnetic field vector in these internal 
layers points increasingly away from the transverse 
direction and this explains the decrease apparent 
in the dashed curve in experiments A and B. On 
the other hand, for experiment C and D the re- 
connection affects a much shallower region of the 
rising tube, and, as a result, the x-component of 
the field as the interface does not decrease. 

Finally, it is found that pressure at the end of 
the experiments, t = 120, has an almost linear de- 
pendence on the value of </>, providing the highest 
magnetic pressure for the cases that do not show 
effective reconnection in Fig. HI The effects and 
importance of reconnection is further discussed be- 
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low. 

4. Magnetic Connectivity 

In the previous section we showed that the 
emerging flux tube reconnects with the coronal 
magnetic field, at least for the cases A - C. As 
a result of this the magnetic pressure below the 
current sheet changes with the relative horizontal 
angle 0o- Thus, it seems plausible that the emer- 
gence process will be strongly influenced by the 
orientation of the coronal field. For example, when 
the two fields are approximately parallel one might 
expect the reconnection process to be slowed sub- 
stantially and the emergence process possibly hin- 
dered. This section investigates the efficiency of 
reconnection, from a global point of view, by 

• measuring the height of the apex and axis of 
the emerging tube in time (Section B~Tj) . 

• measuring the amount of horizontal and nor- 
mal flux that emerges into the corona (Sec- 
tion 



• studying the changes in field line connectiv- 
ity (Section 14. 3p . and finally by 

• measuring the fraction of the tube flux that 
reconnects in time (Section 14. 4p . 

4.1. Height-time relation: apex and axis 

To measure the dependence of the rising mo- 
tion of the tube on the orientation of the ambient 
field, we find the height of the apex and the axis 
of the rising tube as a function of time for the 
experiments listed in Table [U More precisely, a 
large number of field lines are traced from start- 
ing points along the central vertical line. Then 
we find those fieldlines that stay in the tube and 
those that belong to the ambient magnetic field. 
The summit point is then the first point along the 
central line at which the connectivity changes. 

Fig. [5] (left panel) shows the height-time re- 
lation of the apex and the center of the emerg- 
ing tube. The emergence starts with a slow rise 
phase while the flux tu be is below the photosphere 
( Murray et. al.1 ( 20061 ) investigate how this initial 
emergence phases depend on the tube parame- 
ters). This is followed by a rapid rise phase be- 
tween t = 55 and t = 80 during which the apex 




Fig. 5. — Left: The height of the apex (upper 
curves) and the center (lower curves) of the tube as 
a function of time. Right: The associated velocity 
of the apex of the rising system. For line styles see 
Table[U 

of the tube rises through the transition region and 
into the corona. After t = 80 the rise rate slows 
down and settles down to a lower rate that fluc- 
tuates with time and between the various experi- 
ments, but on average follows the same trend. The 
change in the height-time relation around t = 75 
corresponds to the time at which magnetic pres- 
sure balance across the current sheet is achieved, 
as seen in Section 13.11 At this time the width of 
the current sheet shrinks to the numerical resolu- 
tion limit. 

The lower curves in the top panel of Figure [5] 
show that the axis of the tube reaches the photo- 
sphere at t ~ 45 and remains close to the lower 
region of the photospheric layer until the end of 
the simulation. A similar result is fo und in pre- 



vious experiments of flux emergence dFan 




Archontis et al 



2001 



2004; 



Ma gara fc Longcopd 1200 
Manchester et al.ll200 

Finally, the right panel of Fig. [5] shows the rise 
velocity of the apex of the tube. It is seen that 
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all experiments follow the same evolution until 
t « 75, after which a general decrease with super- 
imposed oscillations are seen in the rise velocity. 

4.2. Emerging flux 




50 60 70 80 90 100 110 120 

Time 



Fig. 6. — The graphs show the fraction of flux 
below z = 1.2 Mm as a function of time for the 
five experiments mentioned in Table (pp). Almost 
65% of the flux has emerged into the corona at 
t=120 for experiment A. 

Another global measurement of the amount of 
the emerging flux in time can be obtained by cal- 
culating the amount of new horizontal flux passing 
through the vertical midplane (x-z plane located 
at y = 0) and above a height of z = 1.2 Mm. This 
particular reference height is chosen to be above 
the initial flux tube and below the initial coronal 
magnetic field. In this way, the initial coronal flux, 
independent of its orientation, does not contribute 
to the measurement. 

Alternatively, we can measure the amount of 
horizontal flux that remains below the height of 
z = 1.2Mm in time. This amount of flux is defined 
by: 

/L x pi. 2 
/ B y (x,0,z,t)dzdx. (4) 
-L x J-L z 

Fig. [6] shows normalized by its value at t = 0, 
as a function of time for the five experiments. The 
profiles for the five experiments are almost identi- 
cal until t = 90. This result shows that the amount 
of flux that emerges into the upper atmosphere is 
almost independent of the relative horizontal an- 
gle, After t = 90, we find that the larger the 
initial angle, 0o, is between the emerging flux and 



the coronal magnetic field, the larger is the frac- 
tion of the magnetic flux that emerges into the 
coronal regime at a given time. However, the dif- 
ference in emerged flux between experiment A and 
experiment E is found to be small, close to 7%, at 
the end of the experiment, at which time about 
60% of the initial flux has emerged into the outer 
atmosphere. 

Observa tionally Kubo et al.1 (|2QQ3h ; Spadaro et al 



(|2004l ) and lZuccarello et al.l (|2005l ) have measured 
the development of the normal flux represented 
by the emerging region. They find a time depen- 
dent growth, that over the first few days is close 
to linear. This is followed by a saturation of flux 
and eventually a decrease. The grow rates of the 
magnetic flux in the three cases are naturally dif- 
ferent. One striking difference with the numerical 
experiments is the time scale involved. In the 
observations the timescale of the evolution of the 
system is measured in days, while the experiment 
here only covers about 20 minutes. Therefore it 
will not be possible to make a real comparison, 
but it is interesting to compare the structure of 
the comparable numerical measurement of the 
emerged magnetic flux. 

If one assumes the emerging region is at the disk 
center, then the observed quantity is equivalent 
to integrating the positive (/negative) flux repre- 
sented by the B z component in our experiments: 

= / f V B z (x,y,z,t)dydx. (5) 

Fig. [7l shows the time evolution of the positive 
vertical flux through the transition region (z = 
1.2Mm). The plot shows two different phases in 
the evolution: first a near linear increase of the 
vertical flux until about t = 75, followed by an- 
other phase with a continuously decreasing rate. 
The difference in the time evolution of the verti- 
cal flux between the different experiments is very 
small. 

Despite significant differences in the timescales 
between observations and these experiments, the 
basic structure of the emerging flux appears very 
similar. 

4.3. Field line connectivity 

The magnetic energy released in reconnection 
events is bound to spread along the reconnected 
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field lines due to anisotropic heat conduction. The 
bright structures found in EUV and X-ray obser- 
vations in association with flux emergence, and 
other energy release events, are therefore provid- 
ing vital information about the field line connec- 
tivity in the corona. Comparing the coronal struc- 
tures with models provides us with a possibility of 
understanding the field line structure of dynamical 
events. 

There are several ways to show, in a qualitative 
manner, how the fieldline connectivity changes in 
time for the various experiments. We choose to 
trace field lines starting from one end of the sub- 
merged tube and see if they either connect to the 
other end of the tube or to the corona. A disk is 
selected at the tube end, centered on the initial 
tube axis, and the destination of a large number 
of field lines is determined. Field lines going from 
the one end of the tube to the other are colored 
grey and field lines connecting to the corona are 
colored black. This method indicates the global 
connectivity of the fieldlines between the two flux 
systems. Details of the method can be found in 
IParnell et all (|2004 ). 

Figure [8] consists of five columns. Each column 
corresponds to a different experiment and shows 
how the field lines, which have been traced from 
inside the selected disk, change their connectiv- 
ity in time. The left column, for example, shows 
the connectivity for experiment A. The four pan- 
els in the left column show that the connectivity 
changes first at the outer layer of the disk and 
then moves toward the center following a swirling 



motion. This is because the fieldlines which are 
traced from the outer periphery of the disk, reach 
higher levels in the atmosphere and reconnect first 
with the ambient field. At the end of the simula- 
tion only few fieldlines, which are found in a short 
distance around the center of the disk, have not 
changed connectivity yet. All the other fieldlines 
have already been reconnected with the ambient 
field. 

We also find that at each time the number 
of fieldlines that do not change connectivity is 
smaller, as the relative horizontal angle between 
the two flux systems increases from (/>o = (right 
column) to </>o = 180 (left column). As a result, 
there are very few fieldlines that have been recon- 
nected at the end of the simulation for experiment 
E, whereas most of the fieldlines inside the disk 
have change their connectivity in experiment A. 

It has to be noticed that the disk used for trac- 
ing the field lines is not changed in time, implying 
that the starting points do not represent exactly 
the same footpoints in time. Investigating the drift 
velocity of the flux pattern it is found to be far to 
small to account for the change in connectivity be- 
tween t = 60 and t = 80. A closer examination of 
the field line structure shows that field lines at the 
top of the reconnection s heet continuesly chang es 
connectivity (as shown in lArchontis et al.1 ( 20051 )). 
being slowly pushed towards the flanks of the dif- 
fusion sheet. Here they re close and again become 
part of the flux tube. In other words, flux from 
footpoints at side of the emerging flux region is 
involved in reconnection processes more than ones 
already in the early phase of the emerging process. 
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Fig. 7. — Time evolution of the normal flux of the 
emergence through the z — 1.2 Mm plane for the 
five experiments. For line styles see Table (pQ). 



4.4. Flux connectivity 

To further illustrate the results obtained in the 
previous section, we calculate the amount of flux 
that remains in the tube and does not reconnect, 
normalized to the total flux within the disk, as a 
function of time. 



/black area B y (x,l,z, t)dxdz 
S ( \\ s \ i B y {x,l,z,t)dxdz 



(6) 



<&(t) is shown in the top panel of Fig. [9] for the five 
experiments in Table (pQ). 

Flux emergence through the photosphere starts 
at around t = 60 and reconnection between the 
two flux systems starts shortly after this for most 
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Fig. 8. — Connectivity plots for the five experiments. The five columns represent the A, B, C, D and E case 
respectively. The rows show the connectivity at times close to 60, 80, 100 and 120. The disks show the area 
in the tube at the y— boundary from where the fieldlines are initially traced. 



of the experiments. A measure of the reconnection 
rate is given by the absolute value of the slope of 
the lines in the left frame of Fig. [9j This shows 
that there is a short initial phase where the re- 
connection rate builds up, followed by a period 
of time where the reconnection proceeds with dif- 
ferent, but almost constant, rates in all experi- 
ments. After t ~ 100 the reconnection rate de- 
creases for the three fastest reconnecting experi- 
ments, to a lower level that is roughly maintained 
until the end of the experiments. The reconnec- 
tion rates are simply given by the gradient of the 
connected flux fraction. As already stated ear- 
lier, an unknown amount of flux reconnects more 
than once, implying that gl obal-double-separato r 
bifurcations may take place, Havnes et al. ( 20061 ). 



through which recyc ling of the flux may occur, 



Paranell et al. (l2006h . The estimates of the 



re- 



connection rate are, therefore, only providing a 
minimum value and cannot be used, to estimate 
the reconnection speed in the reconnection pro- 
cesses. This implies that a qua ntitative compari- 
son with lLongcope et"al1 (2005) results is not pos- 



sible. lLongcope et al.1 (|2005h used TRACE and 
MDI observations to estimate the energy transfer 
between a new emerging region and old coronal 
magnetic field. Their results showed that recon- 
nection was not active for a long initial phase, after 
which a large fraction of the emerging flux con- 
nected to the coronal field over a relative short 
timescale. In our simulations, we find a more 
smooth increase in the flux interaction, with a 
clear leveling off towards the end of the experi- 
ment. Hidden in this may well be a significant 
restructuring of the field that is not apparent due 
to limitations in following the connectivity of in- 
dividual flux concentrations. 

Despite this, the graphs still provide informa- 
tion about the general development of the exper- 
iments. The right frame of Fig. [9] indicates that 
experiment A has the highest rate of reconnection 
and that the reconnection rate decreases as the 
angle between the two flux systems becomes less 
favorable and the magnitude of the reconnecting 
field component decreases. This fact is unlikely 
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Fig. 9. — The left panel shows how the flux con- 
nectivity, as defined in Eq. (|6j), changes with time 
for the different experiments. The right panel rep- 
resents the reconnection rate of the experiments 
by simply estimating the gradient of the curves in 
the left panel. The line style is given in Table [TJ 



to be changed since the multiple reconnection de- 
pends on field lines already having changed con- 
nectivity once. Thus, the amount of reconnected 
flux in experiment A increases up to 65% by the 
end of the simulation while it remains close to zero 
when the two flux systems are parallel. Finally, 
we find that the amount of flux that remains in 
the tube at the end of the experiments (t=120) 
decreases nearly as 1 — 0.65 sin(</>o/2). In other 
words, it scales with the orientation, and there- 
fore the strength, of the current sheet. 
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Fig. 10. — Magnetic pressure (thin lines) and gas 
pressure (thick lines) distribution along height for 
the experiments A and D along the central vertical 
line. Line style are given in Table (pQ) 



5. Dynamics of emergence 

In the previous section we showed that the 
height-time relation of the apex of the tube is 
similar in all experiments. It is also found that 
the amount of flux that emerges through a certain 
height as a function of time does not depend on 
the orientation of the ambient field. However, the 



amount of flux that changes connectivity between 
the emerging flux and the coronal magnetic field, 
depends critically on the relative orientation of the 
two flux systems. 

In the experiments where the rate of change of 
connectivity is high, the overlying coronal mag- 
netic field is constantly removed by the reconnec- 
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tion process. Thus, the volume above the rising 
tube is opened and the buoyant system can make 
its way up into the upper atmosphere. On the 
other hand, in the experiments with very little re- 
connection, the coronal magnetic field is not easily 
removed but instead is pushed upwards and keeps 
up resistance to the rising motion of the tube (the 
fieldline topology at the top part of the emerging 
tube is illustrated in Section [6j). Thus, one may 
ask why the rising motion of the tube is not influ- 
enced by the change of connectivity in the different 
experiments. 

In fact, the height-time relation of the apex of 
the tube indicates that the process of flux emer- 
gence is predominantly governed by the dynam- 
ics of the rising magnetized plasma. Thus, in the 
following sections, first we consider the gas and 
magnetic pressure distribution along height, in- 
side the expanding rising volume and across the 
current sheet, for two experiments with different 
initial relative angle (Section 15. then we study 
the temporal evolution of forces that act on the 
upper part of the buoyant tube (Section HT2|) . 

5.1. Pressure distribution 



It has been shown, (Fig. 5 in lArchontis et al 
(2005)), that the three-dimensional current sheet 
which is formed between the two flux systems in 
experiment A, is the location of a rapid change 
in the direction of the magnetic field. In the early 
phase of the evolution of the system the total mag- 
netic field vector goes through a tangential discon- 
tinuity across the current sheet with a clear mini- 
mum at the center of the sheet. As time proceeds 
the direction of the field changes smoothly across 
the interface of the two flux systems following a 
rotational-like discontinuity. 

Figure \W\ shows the gas pressure and mag- 
netic pressure along the central vertical line for 
the experiments A and D at t = 60, 70 and 100. 
The first panel of Fig. [10] (t = 60) shows that in 
both experiments the magnetic pressure is higher 
than the gas pressure inside the expanding vol- 
ume (5 < z < 10) by almost two orders of magni- 
tude. The plasma (3 in this region is therefore very 
low. The magnetic pressure decreases across the 
interface between the two flux systems and has a 
minimum value inside the current sheet. This is 
most easily seen in the top left panel of Fig. [10] for 
the triple dotted dashed line. This position corre- 



sponds to the pronounced minimum of the mag- 
netic pressure occurring at the position of maxi- 
mum electric current, that is due to the tangential 
discontinuity across the sheet. At the same time 
the total pressure has a smooth change over the 
current sheet. This implies that the plasma (3 in- 
creases in the current sheet and becomes larger 
than unity. At this early stage of the experiment 
reconnection at the top of the rising tube has not 
started yet and thus the pressure distribution is 
almost identical in the two experiments. 

The top right panel of Fig. \W\ shows the pres- 
sure distribution when reconnection occurs be- 
tween the tube and the ambient field. In experi- 
ment A, the magnetic pressure still goes through 
a sharp minimum, although the value of the min- 
imum is higher than at t=60. In experiment D, 
instead, then magnetic pressure has a smooth dis- 
tribution across the interface, with no minimum. 
In either case, the gas pressure supplements the 
magnetic pressure across the interface, so that the 
total pressure distribution has no extrema here. 
The reason for the different behavior of the mag- 
netic pressure profiles for the two experiments 
is the presence of an important non-zero, non- 
reconnecting field component in the current sheet 
in experiment D which is nearly absent in experi- 
ment A. 

Finally, the bottom panel shows the distribu- 
tion at a later time, at t = 100. Now, the direc- 
tion of the total magnetic field vector in experi- 
ment A changes smoothly across the current sheet 
following a highly compressed rotational discon- 
tinuity. Thus, there is a finite non-reconnecting 
magnetic component in the current sheet and the 
magnetic pressure there becomes higher compared 
to the magnetic pressure at t = 70. The gas pres- 
sure in the top of the rising system has decreased 
because the dense plasma which was carried up- 
wards has been reconnected. In experiment D, 
the gas pressure and the magnetic pressure do not 
change dramatically compared to the same distri- 
bution at earlier times. Thus, the magnetic pres- 
sure inside the expanding volume in experiment 
D now is much higher than the magnetic pressure 
in experiment A. At all times, the magnetic pres- 
sure exceeds the gas pressure and dense material 
is lifted up against gravity. 

By looking at the pressure profiles with height 
in the bottom panel of Fig. [I0l one notice that 
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in experiments A the current sheet region experi- 
ences a pressure minimum. This is a consequence 
of the Bernoulli effect associated with the lateral 
emission of reconnection outflows from the diffu- 
sion region: the gas pressure depletion in the re- 
gion helps accelerating the plasma from above and 
below into the the current sheet, thereby bringing 
new magnetic flux in to the reconnection site. 

From these plots it is also found that the to- 
tal pressure increases with time in experiment D. 
This is a pileup effect as the upper boundary is 
closed and, thus, the growing excess pressure can- 
not propagate through the top boundary. 

5.2. Forces 

In this section, we focus our attention on the 
forces that act on the upper part of the expanding 
rising volume, just below the current sheet. In the 
following, we consider experiment A. 
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Fig. 11. — Temporal evolution of the vertical 
forces acting below the current sheet along the 
{x = 0, y = 0) line. Shown are the mag- 
netic pressure gradient (dashed-rectangles), the 
tension force (dotted-triangles), the gravitational 
force (dotted-dashed), the Lorentz force (solid), 
the gas pressure gradient (dashed) and the total 
force (thick solid). 

The vertical component of the forces is shown as 
a function of time in Fig. [TU The plot shows that 
the emergence of the buoyant tube is driven by the 
magnetic pressure force that exceeds all the other 
forces at the early stage of the evolution. What is 
also noticeable is that the magnetic pressure force 
and the magnetic tension are equal in magnitude 
and opposite in sign and that they are substan- 



tially larger than either the pressure gradient or 
the gravity force. 
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Fig. 12. — The variation of the magnetic pressure 
force (thin lines) and tension force (thick lines) 
with height for experiment A (triple dotted dashed 
lines) and D (dashed lines) at t = 100. 

The total force is clearly positive accelerating 
the emerging field against gravity until t « 72. 
This is followed by a short period of deceleration. 
After t = 95, all the forces become very small and 
are essentially in balance. This result is consistent 
with the motion of the plasma, which rises with 
an almost constant velocity after t = 95 (see right 
panel in Fig. [5]). Notice that the temporal evolu- 
tion of the total force corresponds well with the 
motion of the apex of the tube shown in Fig. \5\ 
Also, the time (t = 95) at which the forces take 
on very small values corresponds to the time at 
which pressure balance is achieved across the cur- 
rent sheet. 

An analysis of the forces for the experiments 
B-D show that the temporal evolution of the total 
acceleration is similar to the experiment A. This 
result explains why the apex of the tube reaches 
almost the same height at the same time in all 
experiments. To further illustrate the point made 
above we examine the total force for the experi- 
ments A and D for t = 100. The z— component of 
the total force is 

F z = F lz +F g + F p , (7) 

where F g = —pg is the gravitational force, F p = 
— §7 is the gas pressure force and F\ z is the verti- 
cal component of the Lorentz force. The latter is 
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written as 



where the first term describes the magnetic tension 
and gives a force when the fieldlines are curved, 
while the second term is the magnetic pressure 
force that acts from regions of high magnetic pres- 
sure to regions of low magnetic pressure. 




Fig. 13. — Total force as a function of height for 
experiment A (tripe dotted dashed line) and D 
(dashed) at t = 100. 

Fig. H21 shows the distribution of the two terms 
in Eq. (|HJ) along height for the experiments A and 
D for t = 100. On the one hand, the vertical 
component of the magnetic pressure force below 
the current sheet (z < 40) is larger for experiment 
D. On the other hand, the tension force, which is 
a downward force, is also larger for experiment D 
and, thus, the Lorentz force has comparable size 
in both experiments. 

Finally, Fig. [13] shows the distribution of the 
total force in Eq. Q along height for the experi- 
ments A and D for t = 100. The total force at the 
top part of the rising flux system is very small and 
almost identical for A and D. Thus, it seems that 
the total acceleration that acts on the expanding 
rising volume below the current sheet does not de- 
pend on the structure of the overlying field and as 
a result the crest of the tube reaches almost the 
same height at the same time for experiments with 
different orientation of the coronal field. 



6. Topology of the emerging region and 
magnetic reconnection 

The topology of the interacting magnetic field 
is important for understanding the structure of the 
emerging high velocity jets and associated hotter 
plasma distribution. In this section, we show how 
different is the structure of the magnetic field that 
appears in the corona when the relative orienta- 
tion of the interacting magnetic fields changes be- 
tween the experiments. 

The three-dimensional geometry of the cur- 
rent sheet and the jets, emanating from the rims 
of the current sheet, for experiment A an d B, 



have been studied by lArchontis et al.l (|2005l ) and 



Galsgaard et alj (|2005l ). Here, we illustrate the 
projection on horizontal xy planes of the three- 
dimensional structure of the sheet for experi- 
ments A-E. Fig. [TH shows five panels containing 
colourmaps of the total magnitude of the current 
(\J\) on a horizontal cut at a height of 1.7 Mm 
above the base of the corona at t = 100. The 
arrows in the panels correspond to the projection 
of the velocity field. 

The bright patches in Fig. [TH show the loca- 
tion of the highest values of |J|, and correspond 
to the intersection of the horizontal cut with the 
arch-like current sheet. Effective reconnection oc- 
curs and high- velocity outflows are ejected side- 
ways from these sites. The velocities of the jets 
reach values close to 200 Kmsec -1 in experiment 
A and B. In experiments A-D, the direction of the 
jets is aligned with the direction of the ambient 
magnetic field. In experiment E there are no jets 
because reconnection does not occur actively be- 
tween the two flux systems. The arrows in this 
case illustrate the drain of the plasma from the 
uppermost layers of the emerging plasma ball as 
it rises and pushes the ambient field upward. 

The weaker current structure, forming an en- 
closure between the two locations of strong cur- 
rent, outlines the border region between the mag- 
netic flux totally connected to the emerging flux 
tube (inside) and either the reconnected flux or 
the original coronal flux (outside). From this it is 
clear that the horizontal volume of the emerging 
flux tube decreases significantly as reconnection 
becomes more favorably. It is only the flux con- 
tained inside this volume that eventually can end 
up outlining the structure of the emerging flux re- 
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Fig. 14. — Colourmaps of the total current at t = 100. Superimposed is the velocity field (arrows). The 
panels correspond to the experiments A-E. 



gion, while the reconnected field lines connect the 
two flux concentrations to neighboring flux regions 
(here only the corona). 

To further illustrate how the three-dimensional 
reconnection works at the top of the emerging 
flux system we study the topology of the fieldlines 
across the current sheet at t = 100. Figures [15] 
and [16] consist of ten panels that show the field- 
line topology for the five experiments. Three sets 
of fieldlines have been traced from different start- 
ing positions along height and across the interface 
of the two flux systems close to the center of the 
emerging region. The blue fieldlines are traced 
from just below the current sheet and belong en- 
tirely to the rising tube. The red fieldlines are 
traced from just above the sheet and are ambient 
fieldlines. Finally, the yellow fieldlines are traced 
from inside the diffusion region and connect the 
tube with the coronal field. The current concen- 
tration at the interface is visualized with a trans- 
parent isosurface. 

The general picture of the reconnection shows a 
clear difference from the traditional two-dimensional 
configuration. This is because the magnetic field 
vector across the sheet resembles a rotational dis- 
continuity. On the one hand, the uppermost rising 
fieldlines have an orientation which is not perfectly 
aligned with the x axis, as has been explained in 



Archontis et al.l ((2005). On the other hand, the 
orientation of the ambient fieldlines changes from 
experiment A to experiment E so that the rel- 
ative horizontal angle between the two systems 
increases. The product of the reconnection be- 
tween these two sets of fieldlines is another set of 
fieldlines, the yellow lines, which are ejected side- 
ways from the current sheet and establish links 
between the solar interior and the outer atmo- 
sphere. Panels B2, C2 and D2 show that the 
yellow fieldlines at the top of the current sheet 
have an intermediate orientation between the blue 
and the red fieldlines and that they do not stay 
in a two-dimensional plane but they experience 



full 3D-reconnection. The only experiment where 
the initial relative orientation is not favorable for 
reconnection is experiment E. Panels El and E2 
show that there is no reconnection (and thus, no 
yellow fieldlines in the panels). Instead, the coro- 
nal field is pushed upwards and a bended hill-like 
shape interface is formed between the two fields. 
At the late stages of the evolution of the system 
the ambient field slides down along the sides of 
the hill of the emerging flux and some reconnec- 
tion occurs at low heights. 

Through figures [15] and [16] we also get a con- 
firmation of the dependence of the volume of the 
emerged region on the orientation of the coronal 
field that are shown in Fig. [14] More precisely, 
the blue fieldlines illustrate the geometry of the 
outermost layers of the tube and show that the 
expansion, in the transverse direction to the axis 
of the tube, is larger when reconnection is less ef- 
ficient. Indeed, in the case of experiment A the 
outer fieldlines that suffer a large expansion re- 
connect first and as the time goes on more internal 
layers, which expand less and are also less twisted, 
come into contact and eventually reconnect with 
the ambient field. In contrast to this, experiment 
E shows a large expansion of the loops of the upper 
fieldlines as these have not been reconnected due 
to the small angle 0q. This picture confirms the 
indications regarding the location of the current 
structure seen in Fig. [14] 

Finally, if we focus on the shape of the emerging 
fieldlines we find that the blue lines in experiment 
A represent fieldlines which were initially located 
closer to the main axis of the tube and their ori- 
entation was not far away from the y— axis. When 
these fieldlines emerge (see Panel Al and A2) keep 
an almost flat shape in the middle of their crest. 
On the other hand, the outermost fieldlines in ex- 
periments with less reconnection (see for example 
the panels El and E2) represent fieldlines which 
were initially located at the outskirts of the tube 
and they had larger curvature. As these fieldlines 
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rise and expand, they keep their convex shape 
and, thus, the magnetic tension force at the top 
of the emerging tube becomes larger compared to 
the experiments with more efficient reconnection. 
In fact, this has been also shown in Fig.[12j where 
the magnetic tension is plotted against height for 
the experiments A and D at t — 100. 

7. Discussion 

The present experiments were terminated at 
around t=120 due the use of periodic boundary 
conditions in the horizontal direction. For times 
after t=120 the effect of these conditions become 
very apparent on the experiments as the reconnec- 
tion jets have propagated across the system, and 
the subsequent evolution is not showing a freely 
expanding magnetic flux concentration. In the 
Sun, this scenario is possibly more realistic, than 
having open boundary conditions. On the other 
hand, it may cause numerical problems and in- 
fluences the dynamical evolution of the system. 
Thus the timescale of each numerical experiment 
depends also on the boundary effects from the pe- 
riodic conditions. 

Some observations seem to favor a situation 
where a flux rope is emerging into the corona as 
coherent structure dLites et al.lll995h. while others 
([Strouse fe Zwaanl Il999l : IPariat et al1l2004 ) indi- 
cate a pattern where undulating field lines make 
different dynamical effects to release the dense ma- 
terial in the lower parts along them. In our sim- 
ulations, including our previous 3D numerical ex- 
periments, the initial flux tube becomes unstable 
to further expansion into the corona, by building 
up a dominating magnetic pressure force, due t o 
the buoyancy of the tube ( Archontis et al.1 12004 ). 
When the force becomes strong enough it "blows" 
the layers close to the transition region up and into 
the corona, where they rapidly push the overlay- 
ing material away and created a dense, cold mag- 
netically dominated plasma dome. These experi- 
ments show that almost 65% of the normal flux of 
the initial flux tube has emerged into the corona. 
We also find that the axis of the initial tube has 
not fully emerged yet above the photosphere. It 
is wortwhile mentioning that in one of the ex- 
periments we find the formati on of a horizontal 
curren t sheet, first reported bv iManchester et al. 
(|2004 ). which drives internal reconnection of the 



flux belonging to the emerging tube, allowing the 
lower parts of the emerging flux tube to discon- 
nect from the emerging flux system. This provides 
a mechanism to decouple the dense photospheric 
plasma from the field li nes that expand into the 
cor ona (as indicated b y IStrouse fe Zwaanl (|l999l ) 
and lPariat et al.1 ( 2004 )) and provides a possibility 
for forming a structure that looks like an emerged 
twisted flux tube without emerging the entire mag- 
netic flux system. 

As it is seen from the Sections above, the rel- 
ative orientation between the emerging magnetic 
flux and the coronal magnetic field is of great im- 
portance when it comes to the dynamical evolution 
of the flux interaction. The emergence process, 
eventually, produces high velocity plasma jets with 
temperatures in excess of average coronal values 
only in the cases of efficient reconnection. In less 
favorable situations the plasma will not be heated 
much and the spectacular display often seen in 
TRACE movies will not take place. Thus, the 
emergence of an easily observed flux region into 
the hostile coronal environment depends on a rel- 
ative narrow span of angles between the two sys- 
tems. On the Sun, this regime may be increased 
compared to the simple model presented here, due 
to the much larger structural complexity of the so- 
lar environment. 

The experiments discussed here relay strongly 
to the evolution of driven magnetic reconnection. 
Fro m analytical inves t igatio ns of 2D reconnection 
(see IPriest fc Forbesl (|2000f ) and references there 
in), it is found that steady state reconnection de- 
pends critically on the structure of the magnetic 
field, the velocity flow and the value of the mag- 
netic resistivity. This makes it natural to expect 
similar dependences to be carry over to steady 
state 3D reconnection. In relation to this, it is 
obvious that 3D numerical experiments are not 
able to resolve the Reynolds numbers present in 
the coronal plasma. Why should it then be ex- 
pected that the results from the above mentioned 
experiments represent an evolution that may take 
place in the solar corona? Using hyperdiffusion the 
smallest length scale (the thickness) of the cur- 
rent sheet is alway only resolved by a few grid 
points, and it is only at these length scales that 
diffusion of the magnetic field becomes impor- 
tant. Increasing the numerical resolution, implies 
a decrease of the thickness of the current sheet 
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and through this a slight delay in the time for 
the initiation of reconnection. In the driven re- 
connection scenario presented here, the largescale 
of the magnetic field and velocity flow does not 
change significantly when the numerical resolution 
is increased and, thus, the magnetic flux advected 
into the current sheet will remain approximately 
the same. This indicates that magnetic recon- 
nection provides the same new classes of connec- 
tivity. What may change is the local structure 
of the process, where increased numerical reso- 
lution can allow for reaching local turbulence in 
the current sheet as it may become tearing un- 
stable before starting reconnecting. Such a pro- 
cess naturally makes it possible to reach a more 
complicated field line connectivity locally, but it 
also implies that the reconnection process becomes 
independent of the magnetic resistivity and re- 
sponds directly to the amount of flux advec t ed int o 
the current she et (Galsga ard fc Nordlundl Jl996l) ; 
iHendrix et al. I (|l996h : lEvink fc Amid ([20061 )). Fi- 
nally, the choise of resolution for the experiments 
is important. Experiments with low resolution 
yield low values for the energy release obtained in 
the simulations and large values for the timescale 
of the evolution of the system. High resolution, 
especially for the volume occupied by the initial 
tube, may cause problems with the entropy distri- 
bution and affect the buoyancy of the rising tube. 

If the assumption of an emerging magnetic loop 
is correct, and the processes described in this pa- 
per and in our previous work are representative of 
the emergence of flux from the solar interior to the 
outer atmosphere of the Sun, then there are some 
simple features that this model may predict. 

• The process of flux emergence is usually ob- 
served using vector magnetograms, which 
show the appearance of bipolar structures 
at the photosphere. In our experiments, due 
to the highly twisted magnetic flux tube in 
the subphotospheric layer, the orientation of 
the polarity patterns are, in the early phases, 
perpendicular to the main axis of the emerg- 
ing flux tube (North-South orientation). As 
time progresses the two opposite polarity 
spots moving apart in the direction along the 
main tube axis (East- West orientation). 

• The picture of an emerging bipolar region 
with its strong flux region located along the 



axis of the main tube is only true for the 
photospheric layer during the short time of 
these experiments. Slightly above the pho- 
tosphere the topology changes and the ini- 
tial large scale bipolar region structure go- 
ing across the magnetic loop is maintained 
over the remaining time of the various exper- 
iments. Variations between the different ex- 
periments are seen, but the systematic bipo- 
lar patterns going across the main axis of the 
initial flux tube is maintained for all exper- 
iments even for a height that is only 3400 
Km above the photosphere. 

• Magnetic reconnection between the emerg- 
ing and coronal magnetic fields depends on 
the relative orientation between the two flux 
systems. Significant plasma heating and 
high velocity jets will appear when efficient 
reconnection occurs between the two flux 
systems in a coronal environment. For the 
cases where the relative orientation of the 
two flux systems is close to parallel, a weaker 
reconnection may take place closer to the 
photosphere. 

• The emergence rate seems to depend on the 
structure of the magnetic field below the 
photosphere, which in our experiments has 
higher field strength than the coronal mag- 
netic field. The excess magnetic pressure 
therefore expels the flux upwards and, thus, 
opens the volume for a new bipolar structure 
to appear. 

The intense reconnection in the near anti- 
parallel cases makes it possible for most of the 
initial flux of the tube to reconnect with the coro- 
nal magnetic field. If this process continues for 
long enough, then the initial connectivity of the 
flux tube will be totally disrupted. It is therefore 
possible that by emerging a bipolar region into 
the hostile coronal environment, that the initial 
connectivity between the two opposite flux con- 
centrations will be, at the least partly, disrupted. 
How severe this change in connectivity becomes 
depends on several factors, such as the structure 
of the coronal flux, the amount of the emerging 
flux and the duration of the reconnection process. 
A first estimate of the changes in the connectiv- 
ity may be obtained simply by using potential 
models. On the other hand, it seems that these 
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models is not always the best method to follow 
for the study of the time evolution of magnetic 
structures in actively evolving regions, which is 
why we have not used this analysis here. 

Further up in the atmosphere lLongcope et al 



(|2QQ5h have analyzed an interesting emergence ob- 
servation recorded by TRACE. They adopt the 
minimum current model to estimate how much 
flux has reconnected. Using various assumptions 
they find it is possible to account accurately for the 
changes in connectivity. They also find that the 
reconnection between the emerging region and the 
existing coronal flux happens only after some time. 
As the interaction starts, it reaches a high activity 
level that last for a finite time before rapidly lev- 
eling off again. As most of the involved flux seems 
to be accounted for, then it indicates that the sys- 
tem has to build up a significant amount of stress 
before they start interacting. A different possibil- 
ity is that the lower lying magnetic field structures 
are not favorable for the reconnection process and 
that the onset of the reconnection process does 
not start until some of the aligned flux has been 
pushed away and the relative orientation between 
the two flux systems has changed. It is also a pos- 
sibility that the initial weaker reconnection pro- 
cess takes place in regions where the plasma den- 
sity is still so high that the released dissipation 
can not heat the plasma to coronal temperatures. 
Finally there is the possibility that the initial coro- 
nal field in the emergence region is so weak that it 
does not provide for any significant heating as the 
new flux pushes its way up. 

8. Summary 

The emergence of new magnetic flux into an 
existing coronal magnetic field can evolve in very 
different ways. The coronal reaction on the emer- 
gence depends strongly on the relative orientation 
of the two flux systems. Efficient reconnection 
takes place only in the case where the two flux sys- 
tems have almost antiparallel orientations. Taking 
place in the corona, this provides a number of very 
clear signatures, such as, high velocity jets and 
intense heating. The active reconnection process 
clearly affects the volume occupied by the mag- 
netic flux totally connected with the underlying 
magnetic flux tube. 

The reconnection between the two flux systems 



can slow down mainly because of two reasons; Ei- 
ther because all flux in the emerging flux tube has 
reconnected with the coronal field, and by this to- 
tally disrupted the structure of the emerging flux 
rope. Or because the emerging process stops and 
the two flux systems find a mutual balance with 
each other where no, immediate, reconnection is 
taking place. Thus, the length and final structure 
of the new coronal magnetic field strongly depends 
on the total amount of emerging flux, its ability to 
reconnect and on the amount of the coronal mag- 
netic flux. 

For situations where the two flux systems are 
nearly aligned, it will be possible for the emerging 
field to penetrate far into the coronal magnetic 
field before it will create any significant coronal 
signatures revealing its present. These situations 
will only be easily spotted if simultaneous obser- 
vations are made in photospheric and transition 
region lines. 

In our simulations, we have used a simple en- 
ergy equation that does not take into account op- 
tical thin radiation or anisotropic heat conduction 
in the corona. These effects are important ingre- 
dient to include in the experiments before more 
detailed comparisons with observations can be ob- 
tained. These are effects that will be discussed in 
future publication. 

In a subsequent paper we are going to discuss in 
detail the local implications of the ongoing mag- 
netic reconnection process and its implication on 
the dynamical evolution of the different experi- 
ments. 
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Fig. 15. — 3D visualization of the fieldline topology across the current sheet at t = 100 for the experiments 
A (panels Al, A2), B (panels Bl, B2) and C (panels CI, C2). The left column is a side view and the right 
column is a top view of the same snapshot. The current sheet is visualized as transparent isosurface. The 
arrows show the direction of the magnetic field vector. 
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Fig. 16. — Same as in FigUHJbut for the experiments D and E. 
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